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This  final  report  was  prepared  by  Aeronautical  Research 
Associates  of  Princeton,  Inc.,  Princeton,  New  Jersey,  in 
fulfillment  of  USAF  Contract  Number  AF33 (615) -2144;  BPSN; 
5(6899-6l4300l4-000-600-FD);  Project  and  Task  Numbers:  None. 

This  effort  was  done  under  Laboratory  Director's  Funds.  This 
project  is  part  of  a  program  to  provide  the  theoretical  and 
practical  applied  research  necessary  to  maintain  USAF  structural 
design  criteria  for  aerospace  flight  vehicles  and  to  establish 
optimum  structural  requirements,  flight  limitations  and  design 
philosophies  commensurate  with  vehicle  structural  life. 

The  work  was  administered  under  the  direction  of  the  Air 
Force  Flight  Dynamics  Laboratory,  Research  and  Technology 
Division,  Air  Force  Systems  Command,  Wright -Patterson  Air  Force 
Base,  Ohio,  Mr.  Paul  L.  Hasty  (FDTRj,  Project  Engineer.  This 
investigation  was  administered  under  Project  1367  and  Task 
136702;  all  follow-on  work  will  be  funded  under  this  task. 

The  work  reported  in  this  study  was  conducted  by  Aeronautical 
Research  Associates  of  Princeton,  Inc.  with  Dr.  John  C.  Houbolt 
as  principal  investigator,  and  covers  the  period  from  November 
1964  to  March  1966.  The  manuscript  was  released  by  the  author 
in  March  1966  for  publication  as  an  RTD  Technical  Report. 

The  Contractor's  report  number  is  ARAP  Report  No.  83, 

Volumes  I  and  II. 

This  technical  report  has  been  reviewed  and  is  approved. 
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ABSTRACT 


A  development  Is  given  of  four  possible  procedures  for  design¬ 
ing  aircraft  for  gusts  based  on  power  spectral  techniques.  Each 
Is  of  a  different  level  of  complexity  but  all  are  attractive  at 
the  moment.  In  one,  explicit  consideration  of  atmospheric  turbu¬ 
lence  makeup  is  bypassed  altogether.  Another  makes  U3e  of  simple 
over-all  composite  values  of  r.m.s.  value  of  turbulence  intensity 
and  proportion  of  time  in  turbulence.  A  third  is  more  detailed 
and  considers  mission  aspects  explicitly.  A  fourth  is  a  rather 
simple  comparison  approach. 

Items  which  represent  areas  of  uncertainty  are  brought  out, 
chief  among  these  being  the  turbulence  scale  value,  and  the  prac¬ 
tical  means  for  establishing  cutoff  values  in  the  determination  of 
the  basic  frequency  parameter  Nq  of  the  airplane.  The  need  for 
a  more  unified  and  consistent  data  gathering  program  is  brought 
out,  among  various  agencies  and  countries,  and  from  the  standpoint 
of  operational  as  well  as  research  flights.  Analogous  rigid  body 
treatments  by  spectral  and  discrete-gust  lines  are  made  and  com¬ 
pared  so  that  the  discrete -gust  desigr  levels  could  be  used  as  a 
guide  to  help  establish  input  number  by  a  spectral  approach. 

The  material  of  the  report  forms  the  basis  for  the  next 
phase  of  study  which  has  the  aim  of  developing  one  or  more  of  the 
procedures  into  a  firm  design  language.  Computational  data  from 
existing  gust  critical  aircraft  will  help  guide  number  or  design 
boundary  selection  in  this  task. 
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CHAPTER  1 


INTRODUCTION 

Under  contract  with  the  Air  Force  Flight  Dynamics  Laboratory, 
Aeronautical  Research  Associates  of  Princeton,  Inc.  is  developing 
procedures  for  designing  aircraft  for  gust  encounter  based  on 
power  spectral  techniques.  Some  of  the  basic  theoretical  consider¬ 
ations  that  have  evolved  so  far  in  the  first  phase  of  study  are 
presented  in  Volume  I.  Besides  the  treatment  of  gust  design  from 
a  load  exceedance  point  of  view.  Volume  I  also  includes  some 
preliminary  thoughts  on  a  new  concept  for  treating  the  fatigue 
life  of  the  aircraft  structure  as  well,  using  a  common  notation 
and  procedure . 

In  the  present  companion  volume,  further  results  of  the 
initial  phase  of  the  study  are  given.  In  particular,  a  discussion 
is  given  of  four  possible  design  procedures  that  have  been  devel¬ 
oped,  which  range  in  complexity  and  in  degree  of  sophistication 
from  a  simple  lumped  parameter  approach  not  too  unlike  the  past 
discrete  gust  design  procedure,  but  which  incorporate  the  more 
accurate  and  rational  ingredients  of  the  power  spectral  approach, 
to  a  more  detailed  mission  approach  which  takes  into  account  the 
differ  ng  gust  encounter  effects  of  various  stages  of  flight. 

Some  of  the  other  theoretical  work  that  has  been  developed  is 
presented  also. 

Simultaneously  with  the  present  study,  other  related  work  has 
been  carried  on  In  some  of  the  aircraft  companies,  sponsored  also 
by  Wright  Field  and  under  monitorship  by  Dr.  J.C.  Houbolt,  the 
principal  investigator  of  the  present  investigation.  These 
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studies  were  computational  in  nature  and  were  aimed  at  determining, 
for  existing  gust  critical  airplanes  that  were  designed  by  the 
discrete-gust  technique,  certain  basic  response  parameters  that 
are  significant  in  a  power  spectral  approach.  The  intent  in  the 
next  phase  of  study  is  to  continue  the  developmental  work  of  the 
first  phase,  to  incorporate  in  a  design  procedure  sense  the  numbers 
that  were  found  for  the  specific  airplanes  Investigated,  and  to 
end  up  with  a  definite  sp  :tral  design  procedure  (or  several)  with 
definite  design  input  numbers .  One  of  the  underlying  theses  that 
is  being  followed  in  arriving  at  the  proper  or  appropriate  input 
values  is  that  of  applying  the  procedure  in  retrospect  manner  to 
proven  gust  critical  aircraft  with  the  concept  that  if  they  had 
been  designed  originally  by  the  new  procedure,  safe  designs  would 
have  also  resulted. 
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SECTION  2 


SUMMARY  OP  POSSIBLE  DESIGN  PROCEDURES 


A  summary  of  the  possible  design  procedures  that  have  been 
developed  during  the  past  year  of  study  is  presented  in  this 
section^ 

In  Volume  I,  the  development  of  a  "universal"  curve  for  load 
exceedances  due  to  gust  is  indicated.  This  curve  Is  shown  in 
Figure  1;  it  applies  to  any  response  variable  x  of  concern, 
such  as  stress,  bending  moment,  acceleration,  etc.,  where  x  is 
the  increment  in  response  due  to  gust  encounter.  Other  parameters 
are:  N  is  the  number  of  times  per  second  that  the  response  level 
x  is  crossed  with  positive  slope,  NQ  the  number  of  times  per 
second  that  the  response  crosses  the  zero  value  with  positive 
slope,  P  is  the  total  proportion  of  flight  time  spent  in  turbu¬ 
lence,  and  A  is  the  structural  response  parameter  relating  a  , 

A 

the  r.m.s.  values  of  x,  to  ac,  the  r.m.s.  value  of  all  gusts 
encountered  during  flight,  according  to  the  relation 


°X  -  A°c 


The  practical  aspects  of  determining  A  and  Nq  are  outlined  In 
Section  4. 

For  abscissa  values  greater  than  about  10,  the  curve  of 
Figure  1  is  represented  well  by  the  equation 


—  =  (: 
?Nn  \ 


3.36  Act 


6.7^ 
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The  empirical  constants  3.36  and  6.74  are  tentative,  being  based 
on  limited  data  only;  these  constants  are  being  considered  further 
in  continuing  studies  and  should  be  updated  as  more  and  more  oper¬ 
ational  data  are  cast  in  the  universal  form  suggested  here.  In 
the  lifetime  1  of  flight,  the  total  number  of  times  a. given  level 
x  is  exceeded  (crossed  with  positive  slope  for  x  plus,  or  with 
negative  slope  for  x  minus)  follows  from  Equation  (l)  as 


n  =  PTN, 


6.74 


(2) 


This  equation  is  advanced  here  as  a  master  type  equation  for 
designing  aircraft  for  gust  from  a  load  exceedance  point  of  view. 
For  design,  the  level  x  Is  based  on  chosen  design  values,  such 
as  limit  load  or  ultimate  load,  in  accordance  with  the  design 
philosophy  preferred,  while  n  Is  a  stipulated  number  denoting 
the  number  of  times  that  this  level  Is  allowed  to  be  reached  in 
the  lifetime  of  the  airplane  (on  a  statistical  average  basis).  A 
limit-load  reference  will  be  followed  herein,  such  that 


x=Ax=xT7  -  x,  „ 

L.ii.  1-g 


where  x,  _r  Is  the  load  level  in  1-g  flight.  For  the  Implied 
comparison  of  load  exceedances  with  x  on  an  ultimate  load  basis 
consider 


xTT  =  xIT  T  -  x,  „ 
U  U.L.  1-g 
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and  assume  xTr  T 

U  •  ii 

it  follows  that 


1.5  xT  ,  and  let  x  =  xT  =  rx,  Then 
L.L.  L  1-g 


XU  1.5  r  +  .5 


and  in  turn  from  Equation  (2)  that 


This  relation  yields  the  example  results 


or  -i- 

nL 

xl-g 

RU 

1 

107 

1.5 

60 

2 

44 

Thus,  a  single  encounter  of  ultimate  load  conditions  is  equivalent 
statistically  to  on  the  order  of  50  encounters  of  limit  load.  For 
design,  it  would  appear  then  that  a- choice  of  n  =  5  to  10  for 
allowable  limit  load  exceedances  represents  a  good  margin  of  safety. 

From  Equation  (2),  four  design  procedures  of  varying  complex¬ 
ity  may  be  derived,  as  indicated  in  the  following.  The  possibility 
of  treating  fatigue  design  in  analogous  fashion  also  exists,  and 
this  possibility  is  indicated  in  Volume  I. 

2.1.  DESIGN  BASED  ON  STIPULATED  o  .  P,  AND  T  (Presently  recom- 

v 

mended  as  a  second  step  improvement  for  design) 

Equation  (2)  plots  as  shown  in  Figure  2,  where  the  ordinate 
and  abscissa  are  chosen  as  shown  because  of  the  similarity  with 
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Figure  2.  Master  Design  Chart  for  Load  Exceedance  and  Fatigue  Life 


the  familiar  S-N  plots  as  used  in  fatigue  studies  (see  discussion 
on  fatigue  in  Volume  I),  For  design  in  this  procedure,  a  ,  P, 

T,  and  n  are  simply  stipulated  constant  values,  as  for  example 
0=3  fps,  P  =  ,10*  T  =  30,000  hrs,  n  =  5.  Hence  a  convenient 
designation  number  indicating  the  basic  gust  design  characteris¬ 
tics  of  a  'givei  airplane  can  be  established  by  this  procedure; 
for  example,  using  the  example  numbers  just  quoted,  a  3-3000  designa¬ 
tion  airplane  is  indicated,  where  3  refers  to  a  ,  and  where  3000 
refers  to  the  product  of  PT  and  hence  denotes  the  total  number 
of  hours  that  the  airplane  is  expected  to  be  in  turbulence  in  its 
lifetime.  (The  5  is  suppressed  being  chosen  the  same  for  all 
aircraft.)  This  procedure  is  very  much  akin  to  the  discrete  gusi? 
approach  in  which  a  design  gust  velocity  Ude  of  50  fps,  for 
example,  is  stipulated.  The  procedure  has  considerable  merit  due 
to  its  simplicity,  yet  it  still  preserves  most  of  the  basic  ingre¬ 
dients  of  the  problem. 

2.2.  DESIGN  BASED  ON  MISSION  ac,  P,  AND  T  (Presently  recom¬ 
mended  to  study  the  effect  of  extended  flight  times  in  severe 
turbulence;  more  data  along  the  lines  indicated  are  needed  to 
place  this  procedure  on  a  firm  basis.) 

In  this  procedure,  explicit  consideration  is  given  to  the 
various  turbulence  intensity  levels  encountered  and  to  the  amount 
of  time  spent  in  turbulence  at  various  altitudes. 

Figure  3  shows  some  tentative  initial  variations  of  cc  and 
P  that  are  estimated  from  the  consideration  of  existing  data  and 
analysis.  Note  that  only  single  measures  of  intensity  and  propor¬ 
tion  of  time  are  suggested  here;  that  is,  cc  and  P  are  not 
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Figure  3.  Tentative  ac  and  P  Values. 


broken  up  into  severity  categories  as  has  been  done  in  some  past 
studies.  Most  of  the  data  indicate  the  is  nearly  constant 

with  altitude  and  in  the  neighborhood  of  3  to  3.5  fps;  a  constant 
value  of  3.0  has  been  chosen  here.  For  P  the  choice  is  more  in 
doubt;  some  results  sugge-w  curve  a,  others  curve  b.  (The  circum¬ 
stance  is  quite  odd,  since  of  all  quantities  P  should  be  the 
easiest  to  establish.)  At  the  moment,  curve  a  is  felt  to  be  more 
appropriate,  although  proof  of  this  feeling  is  lacking.  The 
recommendation  follows  that  the  curves  of  Figure  3  should  be  up¬ 
dated  from  the  data  that  exist  and  from  those  being  gathered. 

The  effective  value  for  P  applying  to  chosen  missions,  or 
to  flight  paths  which  are  followed  on  the  average,  are  found  from 
the  P  variation  curve  by  appropriately  weighing  the  times  at 
each  altitude  bracket.  Figure  4  shows  example  results  for  the 
assumed  flight  plan  shown  also  in  the  figure,  using  curve  b  of 
Figure  3.  The  3°°  ft/min  climb  and  descent  rate  is  assumed  to  be 
an  average  value  which  takes  into  account  times  spent  in  holding 
patterns . 

With  values  such  as  are  given  by  Figures  3  and  4,  design 
then  proceeds  as  in  the  first  procedure  by  use  of  Figure  2.  An 
illustrative  example  is  given  to  indicate  the  approach: 

Design  Conditions  Given 
x  =  28,000  psi 

T  =  30,000  hrs  «  10,8  x  i07  sec 
n  =  5 

For  The  Mission 

cc  3  fps 
T  =3  hrs  1 

h  I  20,000  ft  /  P  =  -085  fr°m  Flgure  4 
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From  The  Airplane  Preliminary  Design 

A  =  310  pai/fps 

Kq  «  1.1  CPS 

■These  values  yield 


-  .0332 

ptn0  =10.1  x  io6 


which  plots  on  Figure  2  Just  below  the  n  =  3  line  as  shown  by 

the  example  point;  a  safe  design  is  thus  indicated. 

For  special  purpose  missions,  such  as  low  altitude  flight  in 

a  high  percentage  of  severe  type  turbulence,  the  values  of  o 

and  P  must  be  selected  accordingly.  Figure  5  illustrates  how 

a  might  Increase  when  operating  under  more  severe  turbulence 
c 

conditions.  If,  for  example,  flights  are  involved  with  Pg/P,  = .05, 
0-^=3  and  =  7.5#  then  the  ac  for  design  use  is  found  to 
be  3.72. 


2.3.  DESIGN  BASED  ON  COMPARISON  (Recommended  as  a  useful  adjunct 
for  designs  which  bear  close  resemblance  to  proven  gust 
critical  airplanes) 

Write  Equation  (2)  for  airplane  2,  and  also  for  airplane  1, 
and  divide  the  resulting  equations;  the  result  is 


(ptn0) 
( ptnq ) 


(tic 

\  x 


-  6.74 
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which  plots  as  shown  in  Figure  6.  Suppose  that  airplane  1  is  a 
proven  gust-critical  airplane,  and  let  2  be  the  new  airplane 
design.  The  ratios  of  the  PTNq  and  Aac/x  quantities  are  estab¬ 
lished  and  if  the  resulting  point  plots  on  or  below  the  curve  of 
Figure  6  the  new  airplane  is  judged  safe.  This  approach  i3  auite 
attractive  and  has  the  merit  of  simplicity.  Another  advantage  is 
the  use  of  ratios  which  tends  to  minimize  errors  or  uncertainties 
that  may  be  common  in  corresponding  parameters,  such  as  in  the 
choice  of  lift  curve  slope  irr  the  determination  of  the  A’s. 

2.4.  DESIGN  BASED  ON  Nq  VERSUS  x/A  (Presently  recommended  as 
the  first  step  toward  improved  design  procedures) 

In  this  approach,  a  further  simplification  to  the  first 
approach  discussed,  the  values  of  P,  T,  and  a  are  chosen  to 

C 

be  the  same  for  all  aircraft,  and  n  is  specified  to  be  equal  to 
or  less  than  a  certain  number.  On  this  basis.  Equation  (2)  indi¬ 
cates  the  following  very  simple  functional  relation  for  safe  design 

,T  6 -74 

N0  1  C(f)  (3! 

To  establish  the  constant  C,  one  way  is  simply  to  assume  values 
of  P,  T,  a,,  and  n.  Another  way  is  an  empirical  approach; 
thus,  suppose  values  of  Nq  and  x/A  are  computed  for  a  number 
of  existing  gust  critical  aircraft  which  have  proven  themselves 
airworthy  through  many  years  of  successful  operation.  These  values 
when  plotted  would  lead  to  the  shotgun  type  pattern  shown  in 
Figure  7.  A  border  of  the  type  indicated  by  Equation  (3)  may  then 
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Figure  6 .  Gust  Design  by  Comparison  (Note,  stresswise, 
adding  metal  lowers  A). 


Spectral  boundary 


Versus  x/A  Design  Approach  (Points  and 
boundary  are  illustrative) 


be  drawn  along  the  leftmost  points  as  shown,  with  the  inference 
that  design  points  falling  below  this  line  represent  safe  design. 

The  points  and  line  shown  in  Figure  7  are  illustrative  In 
nature  and  do  not  represent  any  specific  computed  values.  How¬ 
ever,  this  approach  is  the  basis  for  the  associated  calculations 
that  are  being  made  for  the  Air  Force  by  some  of  the  aircraft 
companies.  These  calculated  values  of  x/A  and  Nq  are  to  be 
used  in  helping  to  establish  the  definite  position  cf  a  design 
border. 

It  is  to  be  noted  that  all  considerations  of  P,  T,  a  , 

V 

and  n  are  suppressed  or  eliminated  in  this  procedure.  Actually 
the  procedure  bears  a  close  resemblance  to  the  discrete  gust 
design  approach,  but  with  two  basic  and  noteworthy  improvements. 

The  first  improvement  is  associated  with  the  x/A  axis.  Values 
along  this  axis  are  the  equivalent  counterpart  to  Udg,  the  gust 
velocity  value  used  in  discrete  gust  design.  To  see  this,  consider 


Act 


w 


and  multiply  both  sides  by  tj  so  that  T}crx  becomes  equal  to  x^, 
the  design  value  of  then 


By  contrast,  a 
would  appear 


similar  relation  by  th 


l 


discrete  gust  design  approach 
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The  discrete 


indicating  that  Ud  is  similar  in  concept  to  Udg . 
gust  design  procedure  is  thus  represented  in  Figure  7  by  a  fixed 
vertical  line,  as  indicated.  The  newer  approach  discussed  here 
not  only  allows  a  variation  in  but  establishes  the  value  more 

accurately  since  A  is  determined  by  the  more  refined  spectral 
relations  compared  to  the  establishment  of  Adg,  which  is  based 
primarily  on  the  K_  curve. 

o 

The  second  improvement  is  that  of  the  NQ  axis .  In  the 
discrete  gust  approach,  no  account  is  taken  explicitly  of  the 
number  of  times  gusts  load  levels  are  encountered.  The  Nq  axis 
effectively  takes  this  fact  into  account  explicitly.  Thus,  in 
short,  the  boundary  design  curve  of  Figure  7  reflects  that  a 
tradeoff  between  the  A  and  NQ  value  is  possible  in  design; 
specifically,  for  example,  increasing  NQ  must  be  accompanied  by 
decreasing  A. 

In  Section  5,  two  rigid  body  treatments  are  given,  one  along 
spectral  lines  and  one  from  a  discrete  gust  point  of  view.  A 
comparison  of  results  shows  what  terms  in  a  spectral  approach  are 
equivalent  to  terms  commonly  used  in  the  discrete  gust  approach. 
The  material  of  Section  5  is  to  be  used  in  the  continuing  study  as 
a  possible  further  means  to  help  establish  the  position  of  the 
design  boundary  line  on  Figure  7. 
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SECTION  3 

GUST  SPECTRA  AND  L  VALUES 


The  spectra  advocated  here  is  that  associated  with  isotropic 
turbulence  and  due  to  von  Karman,  specifically 


0(1*2)  = 


T 


I  +  |  (1.339  LO)2 
[l  +  (1.339  Ln)2]^" 


(4) 


where  a  is  the  r.m.s.  value  of  gust  velocity,  and  L  is  the 
w 

turbulence  scale.  A  plot  of  the  equation  is  given  in  Figure  8, 

and  an  associated  plot  of  LO0  is  given  in  Figure  9.  This  latter 

figure  is  significant  in  that  it  shows  the  contribution  of  the 

p 

various  frequency  components  to  the  mean  square  value  a  ;  to 
see  this,  write 


a2  =  J  <p{ui)a(in) 
0 


00 


where  dz  =  d(Lfi)/Ln,  or  z  =  log  Lfi.  Thus  with  a  semi -log  plot 

with  Lft  the  log  scale,  the  product  Lf20  becomes  a  distribution 

2  ' 

curve  indicating  the  distribution  of  a  to  the  various  frequency 
components . 

The  value  of  L  to  be  used  in  design  is  at  present  uncertain, 
with  estimates  ranging  from  500  to  5*000  feet.  A  quantity  that  has 
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Figure  9.  Distribution  of  a  with  Frequency 


much  significance  in  establishing  L  is  the  truncated  value  of 
r.m.s.  gust  velocity  defined  by 


*  00 

o  r* 

o£  =  J  <t>{l£l)&(W) 

(“)i 


(5) 


The  insertion  of  Equation  (4)  into  this  integral  leads  to  the 
results  shown  in  Figure  10.  In  application  to  flight  data,  values 
of  o^,  <j  ,  and  L  must  b.  uch  as  to  satisfy  the  curve  of 
Figure  10.  In  general,  is  established  with  good  accuracy 

but  a  is  questionable  because  of  the  presence  of  uncertain  long 
wavelength  components  in  the  gust  velocity  records.  Thus  since 
a1/ow  is  uncertain,  the  deduced  L  is  also.  Stated  in  other 
terms,  for  a  given  a,,  there  can  be  many  combinations  of  a 

-L  n 

and  L  which  when  substituted  in  Equation  (4)  give  </>  curves 
which  fit  the  flight  data  equally  well.  This  fact  can  be  demon¬ 
strated  in  specific  terms  as  follows.  For  Equation  (4),  the 
result  obtained  from  Equation  (5)  for  (LQ)^  >3  is 


.782  oj 


(6) 


With  this  equation  and  the  fact  that  4>(&)  =  L Equation  (4) 
may  be  rewritten 


i 


With  0-,  arbitrarily  chosen  aa  ,003,  this  equation  yields  the 
results  shown  in  Figure  11.  All  curves  have  unit  area  beyond 
Q  ss  a(XXJ  regardless  of  :  or,  if  <k  is  plotted  against  fl, 

p 

all  curves  would  yield  the  same  truncated  value  a^.  Appreciable 
difference  is  noted  in  the  level  of  the  curves  at  lower  C i,  how¬ 
ever,  depending  on  L.  This  lower  range  is  the  range  where  evalua¬ 
tion  of  flight  teat  data  is  difficult  or  Impossible  at  present. 

The  use  of  <t>  curves  in  the  fora  shown  by  Figure  11  may 
become  useful  and  appropriate.  In  general,  the  value  of  A  is 
found  to  depend  strongly  on  the  value  of  L  used,  see  chapter  on 
A  and  Nn  evaluation.  However,  if  the  tf>  curves  of  Figure  11 
are  used,  then  A-^,  the  counterpart  to  A,  and  a3  defined  by 

°x  -  Aow  -  Al°l 

is  found  to  be  practically  insensitive  tc  L.  Thus,  the  use  of  a 

* 

spectral  curve  with  a  specified  truncated  value  may  prove  to  be  a 
worthwhile  direction  to  go,  since  it  minimizes  the  need  for 
establishing  and  designating  a  design  value  of  L.  Possible 
methods  of  improving  the  ability  to  deduce  L  are  being  developed, 
however,  in  the  continuing  research  study. 
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SECTION  4 


A  AND  N0  DETERMINATION 
The  basic  input-output  spectral  relation  is 

<M®)  =  |h(od)  |20(oo)  (8) 

where  cd  -  VO,  With  this,  it  is  desired  to  establish  A  and 

/ 

Nq  as  defined  by  the  following  relations 

°x  "  A0W  (9) 


A 


1 

2 


(10) 


(ID 


Two  questions  arise;  what  value  of  L  should  be  used  in  the  equa¬ 
tion  for  <t>,  Equation  (4),  and  what  should  cd  be.  The  scale  L 
was  discussed  in  a  brief  way  in  Section  3i  some  additional  observa¬ 
tions  are  made  here.  Calculation  studies  show  that,  in  general,  A 
depends  on  L  (cd  fixed).  Typical  results  are  shown  in  Figure  12. 
As  suggested  in  Section  3>  the  use  of  a  truncated  value  of 

r.m.s.  gust  velocity  might  have  merit.  Thus,  use  Equation  (6)  and 


rewrite  Equation  (9)  as 


The  values  for  the  A  values  shown  in  Figure  12  are  also 

shown  in  the  figure.  It  is  seen  that  A^  is  nearly  independent 
of  L. 

Nq  is  found  to  be  insensitive  to  L  ,  Figure  12.  However, 

Nn  is  found  to  depend  on  cu  in  some  cases.  Typical  variations 
of  A  and  Nq  with  o>c  are  shown  in  Figure  13.  A  is  found 
to  level  off  nicely;  however,  Nq  may  level  off,  case  (a),  or 
in  some  cases  may  continue  to  grow  with  increasing  o>  ,  case  (b). 
For  these  later  cases,  some  means  of  establishing  a  practical 
cutoff  frequency  is  needed.  A  rule  of  thumb  advanced  here  is  to 
choose  o>c  at  the  place  where  A  has  leveled  off  quite  well 
and  to  read  off  Nq  at  this  frequency. 

Other  means  for  establishing  cu  are  being  studied. 

v 
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SECTION  5 

SPECTRAL  DESIGN  LIMITS  BASED  ON  DISCRETE  GUST  DESIGN 


A  rigid  airplane  with  a  single  degree  of  freedom  in  vertical 
motion  is  considered  here.  Gust  design  by  the  discrete  gust 
approach  and  also  by  a  spectral  approach  are  considered.  Then  by 
comparison  of  results,  spectral  design  values  are  deduced  from 
the  discrete  gust  design  values  that  have  been  used  in  the  past. 

The  incremental  acceleration  in  g  units  of  the  rigid  body 
due  to  discrete  gust  encounter  is  given  by 


An 


aPpSYe 

2W 


U,  K 
de  g 


(12) 


where  K  is  the  gust  alleviation  factor  which  depends  on 

O 

|ig  =  2W/apcgS,  the  airplane  mass  parameter  (see  Reference  2). 

From  Reference  3*  it  may  be  shown  that  the  transfer  function 
| H(od)  j  for  vertical  acceleration  is  given  by 


where  p.Q  =  4^.  For  use  here,  the  spectral  Equation  (4)  may  be 
put  in  the  form 


5 
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Note  that  care  has  been  exercised  here  to  put  H  and  0  in  forms 
so  as  to  make  subsequent  results  condense  as  much  as  possible  into 
single  nondlmensional  type  curves  rather  than  into  a  family  of 
curves.  The  spectral  equation  for  An  is  then 

*An(k)  =  |HAr,(k)|20(k) 


An' 


Values  of  A  and  Nq  as  found  from  Equations  (10)  and  (ll)  were 
found  to  be 

1 

(13) 


A  =  S£§v  (c\SB 
2W  \2  V 


N0  "  IT  N01 


(l4) 


where  B  and  Nq^  are  shown  in  Figure  l4,  The  near  invariance 
of  B  and  Nq^  to  L/c  is  to  be  noted. 

With  the  parameter  x/A a  as  discussed  in  Section  2  in  mind, 

V 

let  x  be  An;  then  Equation  (12)  divided  by  Ac  gives 

V 


An 
A  a 


P0V, 


eUde 


pVar 


(f) 


B 


(15) 


or  using  Equation  (6)  (with  c  =  a  ) 

w  c 


An  _  PQ^e^de  .885 

Ao°  =  pV01  (kl)^  B 


(16) 


where,  in  analogy  to  the  discussion  of  given  in  Sections  3 

and  4,  here  refers  to  the  truncated  r.m.s.  value  associated 
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i&SaCX 


with  the  area  under  the  4>(k)  spectrum  beyond  k^.  Prom  Figure  14 
and  the  K  curve  of  Reference  2,  the  curve  shown  in  Figure  15 

o 

follows.  The  close  similarity  of  thi3  curve  with  Figure  7  is 
noted.  Thus,  the  results  shown  in  Figure  15,  based  solely  on 
rigid  body  response  considerations,  but  which  considerations  have 
formed  the  basis  for  gust  design  for  many  years,  may  prove  quite 
helpful  in  establishing  the  design  boundary  curve  of  Figure  7.  As 
an  illustrative  indication  of  possible  limits  on  x/A  and  Nq, 


go  back  to  Equations 

(14)  and  (15).  Suppose 

P  = 

P0<  V  =  V  I 

take  Ude  =  50,  and 

choose  Kg/B  =  .26  and 

Noi  ■ 

=  .021  as 

representative  values;  Equation  (15)  indicates 

the 

following 

2L 

An 

An 

(for  a  =  3) 

c 

A 

A  ac 

100 

60.2 

20.1 

200 

76.1 

25.4 

300 

87.2 

29.1 

400 

95.7 

31.9 

while  Equation  (14)  yields 


2V 

c 


20 

30 

40 

50 


.42 

.63 

.84 

1.05 


The  values  of  An/A  and  Nq  are  noted  to  confirm  the  significant 
design  range  shown  in  Figure  7,  and  the  An/A o  values  are  seen  to 

V 

fall  in  the  upper  or  design  range  of  Figure  1.  The  importance  of 
knowing  realistic  values  of  2L/c  and  a  is  brought  out,  however, 

V 
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alternatively  the  Importance  of  knowing  the  truncated  value  as 

used  in  Equation  (16)  is  shown. 


SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


Most  of  the  conclusions  and  observations  arrived  at  in  the 
course  of  study  were  brought  out  in  the  development  in  each 
chapter.  A  summary  of  some  of  the  more  significant  points  as 
well  as  others  not  specifically  covered  is  given  here. 

The  study  has  resulted  in  four  possible  design  procedures  of 
differing  complexity  and  usefulness,  all  of  which  look  good. 
Appropriate  values  of  input  design  parameters,  a  ,  P,  T,  must 

vs 

now  be  established.  Some  of  these  evaluations  must  await  better 
flight  data,  or  at  least  fli.ght  data  reduced  in  a  more  appropriate 
way.  The  form  and  shape  of  the  universal  curve  needs  verification. 
Convenient  analytical  expressions  for  the  universal  curve  are 
being  studied,  and  consideration  is  being  given  simultaneously  to 
establishing  the  associated  and  significant  frequency  distribution 
curves  for  the  r.m.s.  value  of  gust  velocity  a. 

At  present,  much  of  the  flight  data  are  conflicting  and  much 
cannot  be  used  to  advantage  because  the  analysis  was  not  done  in 
a  proper  or  standardized  way  or  because  the  appropriate  quantities 
are  not  presented.  With  respect  to  establishing  a  universal  load 
exceedance  curve,  some  flight  data  indicate  that  the  curve  is 
concave  upward  (on  a  log-log  plor),  as  used  herein,  while  other 
data  Indicate  the  curve  might  be  straight  or  even  concave  downward. 
The  differences  noted  are  believed  due  primarily  to  differences  In 
the  way  various  investigators  take,  interpret,  process,  and  present 
flight  data.  Thus,  a  more  consistent  data  gathering  and  analysis 
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procedure  Is  highly  desirable.  To  this  end,  consideration  is 
being  given  to  the  way  flight  data  ought  to  be  taken  and  analyzed 
to  yield  results  in  general,  or  "universal"  form,  and  convenient 
for  design  application  use. 

The  appropriate  value  or  values  of  L  and  oQ  still  needs 
to  be  found.  Better  means  for  analyzing  flight  data  in  a  meaning¬ 
ful  way  are  being  studied  with  the  aim  of  establishing  realistic 

values  of  L  and  0  .  Likewise,  the  determination  of  Nn  still 

c  u 

presents  some  problems.  Methods  to  place  this  evaluation  on  a 
firm  basis  are  also  being  pursued. 

Other  related  work,  such  as  showing  how  NQ  is  related  to 
the  Initial  curvature  of  the  autocorrelation  function,  or  deducing 
turbulence  velocities  from  mean  wind  profiles,  or  of  combining 
spectral  shapes  to  alter  the  over-all  form,  etc.,  is  continuing. 
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